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Abstract 

We suggest an NMSSM scenario, motivated by dark matter constraints, that may disguise itself 
as a much simpler mSUGRA scenario at tlie LHC. We show how its non-minimal nature can be 
revealed, and the bino-singlino mass difference measured, by looking for soft leptons. 
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I. INTRODUCTION 



The Next-to-Minimal Supersymmetric Standard Model (NMSSM) provides an elegant 
solution to the fi problem of the MSSM by the addition of a gauge singlet superfield S 

mm 

. The superpotential of the Higgs sector then has the form XS{Hd ■ Hu) + \tiS^- 
When S acquires a vacuum expectation value, this creates an effective /i term, /i = A(S'), 
which is automatically of the right size, i.e. of the order of the electroweak scale. 

The addition of the singlet field leads to a larger particle spectrum than in the MSSM: in 
addition to the MSSM fields, the NMSSM contains two extra neutral (singlet) Higgs fields 
-one scalar and one pseudo-scalar- as well as an extra neutralino, the singlino. Owing to 
these extra states, the phej^omeuology of the NMSSM can be s,gmficantly different from the 
MSSM; see Chapter 4 of [5i] for a recent review and references. In particular, the usual LEP 
limits do not apply to singlet and singlino states. Moreover, the singlino can be the lightest 
supersymmetric particle (LSP) and a cold dark matter candidate. 

In this paper, we investigate the LHC signature of an SPSla Gj-like scenario supplemented 
by a singlino LSP. In such a setup, gluinos and squarks have the 'conventional' SUSY 
cascade decays into the bino-like neutralino, ^ ^ B, which then decays into the singlino 
LSP, x5 ~ ^1 plus a pair of opposite sign same-flavour (OSSF) leptons. The X2 decay 
proceeds dominantly through an off-shell slepton. A dark matter relic density of Vth? ~ 
0.1, compatible with astrophysics measurements, is obtained if the x? and/or ^ annihilate 
through pseudo-scalar exchange in the s-channel. 

One peculiar feature of this scenario, taking into account experimental constraints in 
particular from LEP, is that the mass difference between and X2 turns out to be small; 
it reaches at most ~ 12 GeV, and is often much smaller. The leptons originating from the 
bino decay to the singlino, 

(1) 

hence tend to be soft. With the recent interest in soft leptons at the Tevatron Q], searches 
for these should have appeal also beyond the model presented here. 

In the standard SUSY analysis for LHC events, often requiring pri^^) > 20 GeV, there 
is a risk of missing these leptons and wrongly concluding to have found the MSSM instead 
of the NMSSM, with X2 the erroneous LSP and dark matter candidate. Discovery of the 
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Parameter Mi M3 /ies M^^^ M^^ M^^ M^^ M^^ M^^ M^^ M^^ M^^ 

Value [GeV] 120 240 720 360 195 136 133 544 526 524 496 420 521 

TABLE I: Input parameters in for our SPSla-like scenario. The NMSSM-specific parameters A, 
K, Ax and are given in Table HIH 



additional Higgs states will also be very difficult at the LHC in this scenario.^ The aim of 
this paper is to explore the possibility of detecting the X2 ~^ Xi^~^^~ decay, and measuring the 
bino-singlino mass difference, by looking for soft di-leptons at the LHC. Some preliminary 
work on this scenario was carried out in j^. 

In Section [Til we begin by describing the particular NMSSM scenario under investigation, 
and define five benchmark points typifying the small bino-singlino mass difference. We go 
on to discuss the Monte Carlo simulation of these benchmark points in Section UTTl using the 
fast simulation of a generic LHC detector as a basis for studying the possibility of detecting 
the resulting soft leptons in LHC collisions. In Section [IV] we further discuss the extraction of 
mass constraints on the singlino from the di-lepton invariant mass distribution using shape 
fitting, before we conclude in Section |Vl 



II. THE NMSSM SCENARIO 



We use the NMHDECAY 



[10 



ll| program to compute the NMSSM mass spectrum and 



Higgs branching ratios, and to evaluate the LEP bounds: SPHENO 12l is used to calculate 

n n 

the sparticle branching ratios, and MICROMEGAS |13l . Il4j | for the relic density. The SUSY- 
breaking parameters of our scenario^ are listed in Table [B The main difference from the 
familiar mSUGRA scenario SPSla [6] is that we choose a larger Mi = O.5M2 = 120 GeV, 
leading to bino and wino masses of 115 GeV and 222 GeV, respectively, in order to evade 
LEP bounds when adding the singlino and singlet Higgses. In the original SPSla scenario, 
the bino and wino masses are 96 and 177 GeV. The resulting SUSY spectrum, with the 
exception of the NMSSM-specific masses, is shown in Table HIl 



^ The lightest scalar 5*1 and pseudo-scalar Ai are mostly singlet states with masses of about 80-100 GeV 
and decaying dominantly into bb; decays of heavier Higgs states int_o them occur only with branching 
ratios at the permille level. For a discussion of search strategies, see 
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Particle X2 n eR f2 X3 xl X5 *i Ql,r 9 

Mass [GeV] 115 132 143 201 205 222 365 390 397 550-570 721 

TABLE II: Mass spectrum of our SPSla-like scenario; m,-,, „ = mp, n and m~± ~ m,-o . The 
^ 1 ^J'L,R f^L,R Xi^2 ^3,5 

LSP and Higgs masses depend on the NMSSM-specific parameters and are given in Table IIIII 
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1.49 


2.19 
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-49.0 


105.4 


88 
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89 


0.101 
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B 


1.12 


1.75 


-42.4 


-33.6 


112.1 


75 


226 


100 


0.094 


9 X 10" 


-13 


C 


1.20 


1.90 


-39.2 


-53.1 


113.8 


95 


256 


97 


0.094 


1 X 10" 


-13 


D 


1.47 


2.34 


-39.2 


-68.9 


114.5 


109 


259 


92 


0.112 


4 X 10- 


-14 


E 


1.22 


1.95 


-44.8 


-59.1 


114.8 


101 


219 


96 


0.096 


8 X 10" 


-15 



TABLE III: NMSSM benchmark points used in this study. Masses and other dimensionful quan- 
tities are in [GeV]. The other parameters/the rest of the spectrum are/is given in the previous 
tables. 

To obtain a singlino LSP, we further choose A ~ IQ-^ and k, ~ O.IA. This way Xi ~ 99% S, 
and m^o hardly varies with A and k, (Am^o ~ 0.1 GeV). In addition, the trilinear Higgs 
couplings A\ and are chosen such that m^o + m^o ~ m^^ for at least one combination of 



i,j = 1,2, in order to achieve a dark matter density of 0.094 < ^Ih'^ < 0.135 [l5(]. We thus 



The five points used in this study are summarised in Table lllli Points A-E have Am = 
9.71, 3.05, 1.45, 0.87 and 0.60 GeV, respectively. The SM-like second neutral scalar Higgs, 5*2, 

ri 

has a mass of 115 GeV for all these points, consistent with the LEP limit of 114.4 GeV [16i] . 
By contrast, the lightest neutral scalar Si and the lighter pseudo-scalar Ai are mostly 
singlet states, and can hence be lighter than 114.4 GeV. Concerning the efficient neutralino 
annihilation needed to achieve an acceptable dark matter density, for Point A the dominant 
channel is X2X2 ~^ contributing 88% to the thermally averaged annihilation cross section 
times relative velocity, {av) oc l/iVth'^). For Point B, XiXi? X1X2 X2X2 annihilation to 
bb contribute 10%, 15%, and 50%, respectively. Point C has again dominantly X2X2) "while 
Point D has about 50% X2X2 35% X1X2 annihilation. Finally, for Point E, XiXi? X1X2 




obtain a set of NMSSM parameter points with varying Am = m^o — m^o. 
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FIG. 1: pt distributions for leptons from the decay X2 ~^ Xi^^^ in benchmark points A-E. All 
distributions are normalised to unity over the whole momentum range. 



and X2X2 annihilation to bb contribute 29%, 34%, and 13%, respectively. 

Assuming that the non-MSSM nature of the Higgs sector will be very time and integrated 
luminosity consuming to determine at the LHC, or even be difficult to clarify at all, early 
signatures of this scenario will have to rely on the leptons produced by the two X2 ~^ Xi^^^~ 
decays present in the vast majority of SUSY events. Figure [1] shows the pt distributions of 
these leptons for all five benchmark points.^ Clearly, cuts on lepton transverse momentum 
of even 10 GeV will remove the wast majority of events for points B-E, and hence remove 
the one remaining clue to the non-minimal nature of the scenario. However, one should also 
notice that the distributions have considerable tails beyond the simple mass difference Am, 
due to the boost of the X2- Thus, a reduction in the lepton pT-cut holds the promise of 
giving considerable extra reach in this scenario. 

^ For details of the Monte Carlo simulation used, see Section Hill 
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III. MONTE CARLO ANALYSIS 



We perform a Monte Carlo simulation of the benchmark points described above by 
generating our SUSY signal with PYTHIA 6.413 |l7l and SM background events with 



HERWIG 6.510 



18 



19|, interfaced to ALPGEN 2.13 j20| for production of high jet multi 



plicities and JIMMY 4.31 2lJ for multiple interactions. The generated events are then put 



22| 



through a fast simulation of a generic LHC detector, AcerDET-1.0 

Although PYTHIA does not contain a framework for generating NMSSM events per se, 
it has the capability to handle the NMSSM spectrum and its decays. Since our scenario 
predicts the same dominant cross section as in the MSSM, namely gluino and squark pair- 
production, with negligible interference from the non-minimal sector, we use the built-in 
MSSM machinery for the hard process, and take the conservative approach of generating 
only events with squark and gluino production. For the signal, PYTHIA gives a LO cross 
section of 24 pb, and 240 000 events are generated per benchmark point, corresponding to 
10 fb"^ of data. 

For the SM background we have generated a wide variety of samples that in addition to 
two, possibly soft, OSSF isolated leptons at low invariant mass, could potentially yield the 
hard jets and missing energy expected for SUSY events. These consist of n-jet QCD samples, 
generated with jet-parton matching using ALPGEN, and the production of W, Z, WW, WZ, 
ZZ, bb, and ti with n additional jets {n < 3). In addition to this we have also looked at the 
Drell-Yan production of lepton pairs at low invariant masses (3 GeV < mu < 20 GeV) with 
n additional jets. 

These samples were passed through the AcerDET detector simulation. For the scenario 
we consider AcerDET gives a reasonable description of the response of an LHC detector, with 
the exception of soft objects. Given the importance of soft leptons to our study, we therefore 
modified AcerDET as follows: 

1. the pt threshold for leptons was lowered to 2 GeV; 

2. the lepton momentum resolutions used were parameterised from the results of a full 



simulation of the ATLAS detector, as presented in [23l |: 



.3 



^ For muons wc use the results for combined muon system and inner detector tracks with \ri\ < 1.1. The 
electrons are smeared according to a pseudo-rapidity dependent parametrisation. 
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3. we applied parameterised lepton reconstruction efficiencies extracted from results given 
in 

This simulation then incorporates the most relevant effects for the analysis, such as a sensible 
description of the rapidly deteriorating lepton momentum resolution for electrons at pt < 
20 GeV. Also, the applied reconstruction efficiencies fall off steeply for low px, and are 
different for electrons and muons, adding another degree of realism to the reconstruction of 
an invariant mass distribution using both lepton flavors. 

However, there are some issues regarding detector performance at low px that are not 
modelled by these additions, chiefly the introduction of fake electrons through, e.g., the mis- 
identification of charged pions. In particular, one might worry about the potential of pure 
QCD events to fake our signal because of the huge cross section in conjunction with detector 
effects. To improve on the parameterisations used here one would need a full simulation of 
the detector, or even efficiencies from data. We shall show below that all backgrounds with 
pairs of uncorrelated leptons may in principle be estimated from data, assuming lepton 
universality or some knowledge of the degree of non-universality. Therefore, the purity of 
the reconstructed electron sample is less important than the efficiency for reconstructing the 
sample in the ffist place. 



We carry out our analysis along the lines of the 'standard' di-lepton edge analysis 24. l25| . 



see also 
cuts: 



26 



To isolate the SUSY signal from SM background we apply the following 



• Require at least three jets with pt > 150, 100, 50 GeV. 

• Require missing transverse energy fix > max(100 GeV, 0.2Meff), where the effective 
mass Meff is the sum of the px of the three hardest jets plus missing energy. 

• Require two OSSF leptons with pt > 20, 10 GeV. 

After these cuts the background is small compared to the number of SUSY events. For all 
five benchmark points the resulting di-lepton invariant mass distributions have the expected 



Again we use results from combined muon system and inner detector tracks for the muons, where muons 
down to 1 GeV have been simulated. For electrons we use the efficiency of so-called "tight cuts" , defined 



23|, in busy physics events. 
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FIG. 2: Di-lepton invariant mass distributions for point B (left) and point D (right) with stan- 
dard lepton pT cuts. Shown are OSSF (sohd), OSOF (dashed) and subtracted (with error bars) 
distributions. 



edge structure at ~ 80 GeV from the decay chain 



(2) 



and another, much less visible structure, at ~ 110 GeV, due to the same decay through a 
right handed slepton. These can be treated in the usual manner to extract two relationships 
between the four involved SUSY masses, based on the position of the endpoints. 

For benchmark points A and B there is also an excess of events at low invariant mass 
coming from the decays of X2 singlinos, but for points C, D and E the OSSF lepton 
selection cut has removed all trace of the X2 decay. This is demonstrated in Figure [2], showing 
the di-lepton invariant mass distribution for two of the benchmark points. In scenarios like 
C, D and E, one would therefore risk missing the singlino and taking the X2 be the LSP 
dark matter candidate.^ For a further breakdown of the content of the di-lepton invariant 
mass distribution in such scenarios, see j^. 

It is clear that to increase sensitivity to the disguised NMSSM scenario, one needs to 



^ In fact, our SPSla-like scenario is an optimistic one for soft leptons under the standard cuts, in that there 
are extra leptons at hand from the longer decay chain Q to fulfil the cut requirement. 
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1 2 3 4 5 0.5 1 1.5 2 2.5 3 

m„ [GeV] m„ [GeV] 

FIG. 3: Di-lepton invariant mass distributions for point B (left) and point D (right). Shown are 
OSSF (solid), OSOF (dashed) and subtracted (with error bars) distributions. 

lower the lepton px cuts. Due to the hard requirements on jets and missing energy, the vast 
majority of these events should still pass detector trigger requirements. However, lower cuts 
come with the possibility of large increases in SM backgrounds. Most of this background, 
that from uncorrelated leptons, can in principle be removed by subtracting the corresponding 
opposite sign opposite-flavour (OSOF) distribution, assuming lepton universality. However, 
larger backgrounds will increase the statistical error. In addition, a soft lepton sample is 
more vulnerable to the introduction of non-universality from e.g. pion decays. The result 
of lowering the pr requirement on leptons to 2 GeV, after application of the reconstruction 
efficiency, is shown in the left and right panels of Fig. [3], for benchmark points B and D 
respectively. While there is an increase in backgrounds, the effect on the signal is much 
more significant. For both benchmarks, the decay to the singlino is now visible as a large 
excess at low invariant masses. 

To quantify the potential for discriminating between the disguised NMSSM and 
mSUGRA, we show in Fig. Hlthe significance S/ ^J~B of any excess at low invariant masses as 
a function of Am for all five benchmark points and for both the standard and 2 GeV lepton 
Pt cut. The expected number of events 5, in the absence of any signal, is estimated by 
the OSOF distribution plus the expected number of mSUGRA low invariant mass leptons 
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FIG. 4: Significance vs. bino-singlino mass difference Am. We sliow tlie significance both with 
the standard lepton px cut (squares) and the 2 GeV lepton cut (triangles). 



from a fit to the edge at ~ 80 GeV, continued down to low invariant masses. This means 
that the expected number of events can be determined entirely from data. The number of 
signal events S are the events in excess of this. S/ ^/B is evaluated for low invariant masses, 
taking mu < 10 GeV as an upper limit. The exact significance will naturally depend on 
the interval chosen, but 10 GeV should in any case be conservative. At low significance 
there is, as expected, some fluctuation in the significance due to the random nature of the 
signal generation and background generation. From Fig. H] we find that we should be able 
to observe a significant excess down to Am ^ 0.8 GeV, under the assumptions on lepton 
efficiencies described above. ^ However, it is worth noting that even with the standard lepton 
cuts one should be sensitive to mass differences down to 2 — 3 GeV. 



IV. MASS CONSTRAINTS 



In the standard di-lepton analysis the edges at ~ 80 GeV and ~ 100 GeV are used to 
determine the relationship between the neutralino and slepton squared-mass differences, in 



For such small mass differences we may also begin to see neutral displaced vertexes, cf. 
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our scenario m-o — '"^r cind m? — m?o, and the slepton squared- mass m? . With the 

X3 Il/R I-L/R X2 l-L/R 

addition of further edges from longer decay chains, the individual masses m^o^ ^h/R 
m^o can be constrained, although mass differences are determined much more precisely. For 
the SPSla benchmark point, with similar masses to our scenario, one finds that a precision 
of ~ 4% is achievable on the masses of the neutralinos and sleptons involved, when the 



measurement is systematics dominated 27 1. 



In the same manner we could also attempt to extract information on the singlino by 
determining the position of the edge at low invariant masses, giving access to the mass 
difference m™**^ = m^o — m^o, for the three-body decay. Since the shape of an invariant 
mass distribution a priori contains more information than an endpoint, it could be hoped 
that a fit to the whole distribution would further constrain the SUSY parameters involved, 
e.g. setting the scale of the masses as well as their difference {29! . 30 . 31]. In fact, the full 
matrix element for the X2 three-body decay via a virtual slepton, as calculated in {32], is 
used in PYTHIA. From Eq. (11) of [32] we can see that the invariant mass distribution, in 
addition to the neutralino masses, also depends on the left and right handed slepton masses 
and their widths. 

We perform fits to the di-lepton invariant mass distributions at low invariant masses with 



a Gaussian smearing of the shape given in 32], under the assumption that effectively only 



one slepton contributes.^ The Gaussian smearing is meant to emulate smearing by finite 
detector resolution. The results for benchmark points B and D are shown in Fig. [H In 
subtracting the OSOF distribution before fitting, we have taken into account the effective 
lepton non-universality induced by the difference in electron and muon efficiencies. This is 
done by re-weighting pairs of leptons with the inverse of their combined efficiencies, according 
to the lepton momenta involved. This effectively unfolds the non-universality effects on the 
invariant mass distribution from the differing efficiencies, at the cost of increasing statistical 
errors due to large weights. It should be a simple extension to include geometry dependent 
efficiencies into this re-weighting. Naturally, the re-weighting can only be effective if the 
errors on the measured lepton efficiencies are small compared to the other errors involved in 
the fit. The resulting differences in shape can clearly be seen by comparing Figs. [3] and [51 



^ For our benchmark points the right handed slepton contributes ^ 90% to the decay ampHtude. We have 
checked that the fits are completely insensitive to whether there are one or two sleptons participating. 
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FIG. 5: Di-lepton invariant mass distributions for point B (left) and point D (right) after re- 
weighting with lepton efficiencies. Fits (in red) are described in the text. 



Benchmark point A B C D 

m^o-m^o 9.77 ±0.03 2.98 ± 0.02 1.39 ±0.03 0.92 ± 0.02 

mr-m^o 46.5 ± 12.7 52.7 ±21.9 69.0 ± 53.6 57.2 ± 95.8 

X^/ndf 1.20 1.29 2.06 0.90 

TABLE IV: Mass differences in [GeV] and fit quaUty x^/ndf from fits to the di-lepton invariant 
mass distributions. 



Despite our hopes, we find that the shape fits do not constrain the slepton width at all, 
nor do they constrain the absolute mass scale significantly. For benchmark point A, with 
the largest statistics, it indicates a singlino mass of m^o = 83.2 ± 44.1 GeV. However, when 
parametrised in terms of the bino-singlino and slepton-singlino mass differences, keeping 
also the scale as free parameter, the fits give quite good bounds, which can be found in 
Table IVl 

For both benchmark points A and B, with the larger statistics, the fit gives a useful bound 
on the slepton-singlino mass difference. These can be compared to the nominal values of 
mi — rriyO = 37.1 GeV and 30.4 GeV, for points A and B respectively. This sensitivity can 



12 



be understood physically as the effect the proximity of the slepton pole has on the shape of 
the invariant mass distribution. For much larger slepton masses this sensitivity should go 
away. For all four points the fit gives an accurate determination of Am. Comparing to the 
nominal values given in Section [TTl in particular the result for benchmark point B indicates 
that there are potential sources of significant systematic error, larger than the statistical 
errors with 10 fb~^ of data. 

With the information obtainable from a long decay chain involving an on-shell slepton, 
see above, the absolute singlino mass can be found with the same precision as the other two 
neutralinos involved, i.e. around 4% for our scenario, meaning that we are dominated by 
the errors of the long decay chain. The results on the slepton-singlino mass difference in 
Table HVl would indicate, for benchmarks A and B, that the decay ([T]) occurs dominantly 
through a different slepton than the decay ([2]), which one may speculate could in turn give 
some restriction on the neutralino mixing parameters. 

V. CONCLUSIONS 

We have demonstrated that lowering the requirements on lepton transverse momentum 
in the standard search for the SUSY di-lepton edge may reveal unexpected features, such 
as the NMSSM in disguise. While our numerical results are sensitive to the exact lepton 
efficiencies and momentum resolutions at low transverse momenta — to be measured by the 
LHC experiments — the OSOF subtraction procedure ensures that the background can be 
estimated from data and that the NMSSM scenario in question is both discoverable down to 
very small bino-singlino mass differences. Am = m^o — m^o ~ 0.8 GeV, and that this mass 
difference is measurable to good precision. We have also shown that the di-lepton invariant 
mass distribution has some sensitivity to the slepton-singlino mass difference. 

We would also like to note that, since virtually all SUSY cascades in these scenarios will 
contain two decays of the type ([1]), this lower edge in the di-lepton distribution may appear 
much earlier than the 'standard' decay through a slepton, if at all present, provided that 
the soft leptons are searched for. This may in fact be an early discovery channel for SUSY. 
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